Myoelectric Control of Artificial Limb by Quantum Information Processing 
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Precise and elegant coordination of a prosthesis across many degrees of freedom is highly desired 
for rehabilitation of people with limb deficiency. Processing the electrical neural signals, collected 
from the surface of the remnant muscles of the stump, is a common way to activate certain function 
of the artificial limb. Based on the assumption that there are distinguishable and repeatable signal 
patterns among different types of muscular activation, the problem of the prosthesis control reduces 
to the pattern recognition. Widely accepted classical methods for pattern recognition, however, can 
not provide simultaneous and proportional control of the artificial limb. Here we show that quantum 
information processing of the neural signals allows us to overcome above difficulties suggesting a 
very simple scheme for myoelectric control of artificial limb with advanced functionalities. 

PACS numbers: 87.85.E-, 87.19. Ir, 03.67.Ac 



Information embedded within the electrical neural 
pulses, controlling muscle contractions in our body, can 
be extracted from surface or intramuscular myoelectric 
signals, which are typically summarized into so-called 
electromyogram (EMG). While acquisition of intramus- 
cular signals could cause serious ethical concerns and 
may lead to infection due to its invasive procedure, for 
the last forty years, surface EMG signals have been the 
only source powered prosthesis control [Tj. Currently, 
the commercial prostheses utilize simple processing of 
the surface electromyogram, and can provide very lim- 
ited functionalities To improve the functionality of 
myoelectrically controlled prosthesis, pattern classifica- 
tion algorithms for surface EMG have been extensively 
investigated in the academia. It has been shown, for 
instance, that with properly selected features and clas- 
sifiers, one can achieve very high classification accuracy 
(more than 10 classes of movements with less than 5% of 
classification error) 

The academic success of myoelectric control based on 
pattern classification, however, has not translated into 
significant commercial and clinical impact as one would 
have expected. So far, none of the commercial prosthe- 
ses is using pattern classification based controller. One of 
the main problems with the pattern classification for my- 
oelectric control is that it leads to very unnatural (for the 
user) control scheme While natural movements are 
continuous and require activations of several degrees of 
freedom (DOF) simultaneously and proportionally, clas- 
sical schemes for pattern recognition allow activation of 
only one class that corresponds to a particular action in 
one decision, i.e. sequential control. Moreover, all these 
classes as well as their superpositions must be previously 
learned. Simultaneous activation of two DOFs is thus 
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recognized as a new class of action, but not as a combi- 
nation of known actions. This means higher complexity 
of the classifier and consequently less robustness. In ad- 
dition, this also implies that the user must spend more 
time in training and learning, resulting in higher rehabil- 
itation cost and more frustration. 

Recently there are attempts in addressing the issue 
of simultaneous and proportional control of prosthesis. 
Instead of taking classic pattern recognition approach, 
for example, multilayer perceptron neural network was 
used to estimate joint force ^ and joint angles jS^ of 
the actions from surface EMG. One of the main limi- 
tations of this multilayer perceptron approach is that it 
still requires the data from combined activation of multi- 
ple DOFs in the training set. Alternatively, non-negative 
matrix factorization has been used to estimate the joint 
force |3] and the joint angles [Sj. This approach only 
requires data from single DOF activations for calibra- 
tion, but has inferior performance in comparison to mul- 
tilayer perceptron network approach. It is worth not- 
ing that both mentioned control schemes require solving 
optimization problems during training stage, which in 
turn demand significant computational power for multi- 
ple DOF. The training, moreover, needs to be performed 
on large data sets with typical size of lO"* — 10^ sam- 
ples, which need to be acquired during time-consuming 
sessions. Taking into account that the classical con- 
trol schemes, such as multilayer perceptron, linear per- 
ceptron, linear discriminant analysis, Gaussian mixture 
model and hidden Markov model, do not show significant 
difference in classification power [7], in this paper we de- 
velop a radical approach to the problem of myoelectric 
control based on quantum information processing. 

During the last two decades, information encoding into 
the states of quantum systems and its processing accord- 
ing to the laws of quantum mechanics have been demon- 
strating impressive advantages over classical information 
processing [5]. Typically these advantages are discussed 
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FIG. 1: The three DOF of the wrist to be emulated by the 
prosthesis. 



in the context of computational complexity [5] and com- 
munication security |10| . Recently, however, we have 
shown that quantum information processing can dramat- 
ically increase the capabilities of the simplest learning 
machine - perceptron It has been demonstrated, in 
particular, that quantum perceptron is able to learn arbi- 
trary logical functions, perform classification of data with 
an overlap and even conclude when the given features do 
not belong to previously seen classes - an unattainable 
tasks for its classical counterpart. Moreover, quantum 
perceptron learning rule does not involve any optimiza- 
tion procedure that makes the quantum information pro- 
cessing immune to the curse of dimensionality, which se- 
riously limits power of classical learning models |I2) . 

In this paper we adopt the concept of quantum per- 
ceptron to provide simultaneous and proportional myo- 
electric control of the artificial limb. While superposi- 
tion principle of quantum mechanics allows us to recog- 
nize combinations of actions without their previous learn- 
ing, the fact that quantum state is normalized, i.e. 
{tjj I'ijj) = 1, ensures proportional control of the prosthesis. 

Let us focus on a particular case of limb deficiency - 
transradial amputation, or the amputation at the fore- 
arm. Our choice is mainly based on two factors. First of 
all, this type of limb deficiency represents a large portion 
of the upper-limb amputations. Second, wrist movements 
are complex enough and require activation of multiple 
DOFs. In our analysis we consider a wrist prosthesis 
with three DOF, namely flextion-extension, radial-ulnar 
deviation and pronation-supination, as illustrated in Fig- 
ure [1] 

Let us suppose that n electrodes are placed on the 
surface of the deficient limb. Particular choice of the 
location of the electrodes strongly depends on particu- 
lar amputation and remnant musculature, therefore it 
is discussed elsewhere. There are four classical features 




FIG. 2: The principal scheme of the myoelectric control of 
the artificial limb. The electrodes placed on the surface of 
the deficient or normal (as in the left picture) limb detect 
the electrical neural signals, which are summarized in EMG 
(two particular examples of EMG compound of seven elec- 
trode channels are displayed upward). These figures repro- 
duced from 5 , with permission from Biomed Central. EMG 
features are encoded into the amplitudes of quantum states. 
During the learning stage, the three sets of positive operator 
valued measurements (POVM) operators are deduced, which 
control the three DOF. These operators are utilized to decide 
on kind of action to be performed by the prosthesis during 
the autonomous operation stage. 



that can be extracted from the raw signals of each elec- 
trode channel: mean absolute value, zero crossing, slope 
sign change and wave length [13] . Thus, for each of the 
four features the feature space is of dimension of n. Let 
one of these features, mean absolute value for example, 
be encoded into the states of a (discrete) n-dimensional 
quantum system, so that the component a.i \i) of the state 
lip) = ^ |i) represents the mean absolute value of the 
signal from the ith electrode. Here, the amplitude of 
the signal ai is normalized over amplitudes from all the 
channels by factor l^iP to fulfill the normalization 
condition — 1- Let the quantum state be the 

input of a single-layer network of quantum perceptrons. 
In the most simple architecture that we are going to con- 
sider, the number of the perceptrons in the network is 
equal to the number of DOFs to be controlled, so that 
each quantum perceptron governs just one specific DOF. 
Without loss of generality, let us suppose that the per- 
ceptron Dl controls flextion-extension, D2 ~ radial-ulnar 
deviation, while D3 - pronation-supination. The control 
scheme is represented schematically in Figure [2j 

Before free control of the artificial limb, the subject 
with the limb deficiency is instructed to focus on per- 
forming particular actions to learn how to control the 
prosthesis. Technically, the aim of the learning stage is 
to deduce the association between the EMG features and 
the joint kinematic. For example, given a corresponding 
command, the subject is trying to perform flextion acti- 
vating DOF Dl. EMG features, extracted from received 
myoelectric signals, are encoded into the quantum state 
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I/) = J2fi\^)- Since efficient control of artificial limb 
demands not just recognition between particular classes, 
but also prediction of a degree of chosen action (angle or 
force), the learning must be repeated for different angles 
of flection. The quantum states, constructed from EMG 
features corresponding to different angles 6i are combined 
into POVM operator PjP^ = |E "i/^) (E where 

Following the above procedure, learning extension 
leads to the construction of operator P^^. The operators 
and P^^ summarize all training data that activate 
the Dl. However, these operators do not necessarily form 
a legitimate POVM. The fundamental property of POVM 
operators is that they form a complete set, therefore, we 
need to define the third operator P^^ = / - Pf^ - P^^ 
to fulfill this condition. Here / is the identity operator. 
Operator P^^ thus collects all features that correspond 
to EMG activities, which do not lead to activation of 
the Dl. In the same way as described above, the learn- 
ing procedure is to be repeated for the other two DOFs 
D2 and D3. As the result of the learning stage we have 
deduced nine operators that govern the three DOF ac- 
cording to the classification of the selected feature (mean 
absolute value in our case). The learning must be re- 
peated for the other features the same way as discussed 
above. 

Let us now see how the trained network of quan- 
tum perceptrons responds to the EMG during subject's 
autonomous control of the artificial limb. Desiring to 
perform an action, the subject generates neural signals 
which, as before, are acquired from the limb surface, ex- 
pressed in EMG and encoded into a quantum state j-^)- 
Each quantum perceptron computes the expectation val- 
ues {ipl P lip) for the nine operators |15) . The expec- 
tation values are interpreted as the strength of activa- 
tion. For Dl, for example, the three expectation values 
are computed as f — {ip \ Pj'^ j?/)), e = (i/'l Pe^ IV') and 
^Di — (V'l Pq^ The expectation values / and e con- 
tain information of how likely the given EMG correspond 
to flcxtion or extension. The expectation value O^ia gives 
the probability that the flextion-extension DOF is inac- 
tive. If / > e, then the prostheses makes flextion on the 
degree J^pcfpniv where 0/ is the maximal flextion 

observed in the training data and Tr(...) stands for trace 
operation and gives the overlap between the operators 
P^^ and Pg^^. If, in contrast, e > /, the prosthesis per- 
forms extension on ^ r^^jDipmy where 0e is the maxi- 
mal observed extension. Here, we assumed that a linear 
combination of DOFs corresponds to a linear combina- 
tion of corresponding features. This assumption is valid 
for mean absolute value, although may not be true for 
an arbitrary chosen feature 

As the result of the analysis of the given state IV'), 
which encodes mean absolute value feature, the network 
of three quantum perceptrons returns three angles to the 
mechanical system of the prosthesis. These angles de- 



fine the simultaneous and proportional action completely. 
Interestingly, there is additional control information en- 
coded into the three expectation values 0di,^D2 and O^o. 
As we have noted before, the expectation value Qdi tells 
us the probability that the state \ip) does not represent 
neither flextion nor extension. This may be interpreted 
as the given signal activates the other two DOFs, i.e. 
Odi = D2 -t- D3. Such interpretation gives us three ad- 
ditional equations with three unknowns. Resolving this 
system of linear algebra equations, we can deduce addi- 
tional proportionality between the activation of the three 
DOFs. 

We tested the proposed control scheme on a set of 
EMG acquired from an able-bodied subject who per- 
forms wrist contractions. Eight electrodes were placed 
around the circumference of the forearm, with equal 
inter-electrode distance j4]. It is important to note that 
even for a single DOF there are two types of movements: 
direct action from the rest position and return action to 
the rest position. EMG for these actions may differ signif- 
icantly, although formally correspond to the same spatial 
angles. In our analysis we always use direct actions both 
for the training and recognition of new patterns. 

Our training data set collects only those movements 
of the subject that activate one of the three degrees of 
freedom. The test data set contains complex movements 
without any restrictions. The raw EMG that correspond 
to all these movements were recorded with 1024 Hz sam- 
pling rate. Mean absolute value feature was extracted 
from the EMG averaging the signals over 100 ms win- 
dow. In our analysis we focused on just two DOFs Dl 
and D3, as they are the most desired functions not yet 
have commercial availability from trans-radial amputees. 
The purpose is to test whether our control scheme can 
recognize the combinations of these movements, being 
trained on just single DOF activations, i.e. the system 
can extrapolate automatically from single DOF data to 
their arbitrary combinations. 

To test the control scheme, we initially used small sets 
for training consisting of just 500 data strings for each 
action, i.e. flextion, extension, pronation and supina- 
tion. Each data string consists of the mean absolute val- 
ues acquired from the eight channels and the respective 
joint angles of the DOF activation. The testing data 
contains 8216 data strings divided on 55 blocks; each 
block encodes a multiple DOF activation in a certain 
angle range. In 15 of the 55 blocks, the classification 
error appeared: 5 mistakes in flextion-extension, 9 - in 
pronation-supination and once both DOF were misclas- 
sified. To analyze the accuracy of classification, in cases 
when it was successful, we used so-called performance in- 
dex, which is widely used as a global indicator of quality 
of the estimation [S]. The performance index for an fc"^ 
DOF is given by 



where Ui it the joint angle of the k DOF, is the 
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corresponding estimate, Oi is temporal average of and 
the summation is to be done over all data samples. Sim- 
ilarly, the global performance of the estimator is defined 
through the sum over all K DOFs as 

i?2 = 1 - Sf^i E - (2) 

The performance of flextion-extension recognition is 
found to be 0, 834, while the performance of pronation- 
supination classification is 0.224. It is not surprising that 
the performance for D3 is much lower than for Dl. Pre- 
vious studies also reported lower performance on D3 than 
Dl in H]. This is mainly due to the fact that mus- 
cle responsible for D3 are deep muscles. As a results, 
their EMG are easily masked by the EMG of superfi- 
cial muscles, such as flexor muscles and extensor mus- 
cles. This masking effect is particularly pronounced dur- 
ing combined activations of Dl and D3. Nevertheless, the 
global performance of the two degrees of freedom recogni- 
tion 0, 715 is comparable to the performance of classical 
schemes Illic- 
it is important to check, how the efficiency of the esti- 
mation changes with the growth of the training set size, 
because using more data for the training is a standard 
way to improve performance of the model. We subse- 
quently used large data sets for the training: 2000 sam- 
ples for each action. For testing we used the same 55 
blocks of unseen data. In result, in 10 blocks errors were 
detected: 3 for Dl and 7 for D3. Moreover, the global 
performance increased to 0, 736. Thus, the efficiency of 
classification increased with the size to the training set. 

It is important to stress on the role of the overlap be- 
tween operators P^^ and P^^ and between operators 
PjP^ and P^^ in estimation of the respective DOFs. In- 
deed, these pairs of operators, by their construction, dif- 
fer specifically in those components that are crucial for 
activation of the corresponding actions. It has been ob- 
served that the overlaps vary in certain range with the 
growth of training data approaching (on the large scale) 
some optimal value. As closer the overlap to the opti- 



mal values as higher the performance of estimation. In 
practice, closeness of the overlaps to their optimal values 
would indicate that further learning is not necessary. 

Finally, we would like to admit that there are many 
nontrivial ways to further improve our control scheme. 
For example, it may be beneficial to include more quan- 
tum perceptrons in the network, so that several percep- 
trons control one DOF (direct and return actions, for 
example). It is also possible to construct a database of 
POVM operators out of the training data, so that one 
set of operators controls movements in a particular an- 
gle interval. A particular set of operators can be called 
from the database depending on overall intensity of the 
EMG. Also, different features may be encoded and ana- 
lyzed simultaneously due to the fact that quantum am- 
plitudes are complex numbers in general. Finally, our 
control scheme can be combined with classical control 
schemes, since its implementation does not demand any 
optimization. 

In conclusion, we showed that quantum information 
processing can realize intuitive simultaneous and propor- 
tional prosthetic control over multiple DOFs, with only 
training on individual DOF data. We suggested a simple 
control scheme for the wrist prosthesis having advantages 
of optimization absence during the learning stage and not 
demanding large training data sets (size of which is pro- 
portional to the time spend by a subject in a lab and by 
implication his/her frustration in the procedure). 
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\a) = s-'id 1^) = represented through 

quantum states of a n-dimonsional quantum system and 
we know that the vectors belong to one class. Let us 
define the operator P = j^Qa+bfia+b]) • Being a part of 
a POVM set, this operator ensures nonzero probability 
of outcome for any given linear combination of vectors 
Q |a) + /? \b), where \a\'^ + \I3\'^ = 1. Taking into account 
that statistics of future inputs can not be deduced, it is 
reasonable to accept that all superpositions of the vec- 
tors I a) and \b} can be observed with equal probability. 
The operator P maximizes the average probability of de- 
tection of an arbitrary superposition of the two vectors; 
thus it is optimal POVM clement by construction. This 
consideration generalizes straightforwardly to the case of 
an arbitrary rmmbor of feature vectors. If two feature 
vectors \a) = X)'**!*) ^^'^ W = represent the 

same motion but on differen angles 6a and 6b respec- 
tively, the operator P = \aaa + abb) {uaa + abb\, where 



aa = 6a/ {9a + 6b) and at = 6b/ {6a + 6b), takes into 
account contribution of the feature vectors to the activa- 
tion of the motion. This definition of the operator P is 
not unique. It leads, in particular, to the normalization 
of the all action angles to the maximal angle observed in 
the training data. 
[15] A single measurement with the POVM on a given fea- 
ture vector implies three possible outcomes, namely the 
feature vector belongs to either the first or the second 
class, or the feature vector is essentially different from 
the two classes. However, the result of the single mea- 
surement is of no use for myoelectric control of artificial 
limb, because it would lead us to a choice between the 
two actions and no action. The expectation values of the 
operators, in contrast, are irmch more informative: they 
give us probabilistic interpretation what kind of action 
the given feature vector corresponds to. 



